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Abstract 

Coherent control of femtosecond two-photon absorption in the intermediate-field regime is an- 
alyzed in detail in the powerful frequency domain using an extended 4*'*-order perturbative de- 
scription. The corresponding absorption is coherently induced by the weak- field non-resonant 
two-photon transitions as well as by four-photon transitions involving three absorbed photons and 
one emitted photons. The interferences between these two groups of transitions lead to a difference 
between the intermediate- field and weak- field absorption dynamics. The corresponding interfer- 
ence nature (constructive or destructive) strongly depends on the detuning direction of the pulse 
spectrum from half the two-photon transition frequency. The model system of the study is atomic 
sodium, for which both experimental and theoretical results are obtained. The detailed under- 
standing obtained here serves as a basis for coherent control with rationally-shaped femtosecond 
pulses in a regime of sizable absorption yields. 
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I. INTRODUCTION 



Femtosecond pulses offer unique ways to coherently control photo-induced quantum dy- 
namics of matter PQ [21 El IH [S] ■ The corresponding key characteristic is their coherence over 
a broad spectrum. Multiphoton absorption processes in atoms and molecules, which are of 
fundamental scientific importance as well as applicative importance to the fields of nonlinear 
spectroscopy and microscopy, are among the processes that have been controlled most effec- 
tively by shaped femtosecond pulses [5l[6l[71|8l[9l[l0l[IIl[T2l[T3l[Tl[T5l[T6l 
The control principle is the coherent manipulation of interferences among the manifold of 
initial-to-final state-to-state multiphoton pathways that are induced by the pulse. Con- 
structive interferences lead to absorption enhancement (i.e., increased transition probabil- 
ity), while destructive interferences lead to absorption attenuation (i.e., decreased transition 
probability). The interference manipulation is implemented by shaping the femtosecond 
pulse [21], i.e., manipulating the spectral phase, amplitude, and/or polarization of its dif- 
ferent frequency components. Hence, in order to fully utilize the coherent control potential 
of a given excitation scheme, the ideal line of action is shaping the pulse based on an 
initial identification of the different multiphoton pathways and their interference mecha- 
nism. When such identification is not possible, a practical partial solution is to use auto- 
matic experimental optimization of the pulse shape using learning algorithms that generally 
considering the system as a black box [22] • The lesson, which can be learned from the 
many successful coherent control studies of multiphoton absorption conducted in the past 
[5l|6l[ai8l[3IIOllIIlll2l[T3l[Tl|15lll6], is that this ideal line of action is feasible and 
very powerful once the photo-excitation picture is available in the frequency domain. This 
is possible only within the framework of perturbation theory, where a valid perturbative 
description in the time domain is (Fourier) transformed to the frequency domain. How- 
ever, until recently the frequency domain has been exploited only in the weak-field regime 
[51El[ZlEl[Sl[ini[IIl[I21[I31[Il[I5], where the N-photon absorption is described by per- 
turbation theory of the lowest non- vanishing order, i.e., the N*'* order. Physically it means 
that the N-photon absorption is coherently induced by all the possible initial-to-final state- 
to-state pathways of N absorbed photons. For two-photon absorption the lowest order is 
the 2"^^ one, involving all the pathways of two absorbed photons [3 |6l [71 |8l |9l [TOl [11] . The 
downside of being limited to the weak-field regime is the low absorption yields associated 
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with it. For two-photon absorption they are typically below 0.1% population transfer. 

Recently [T^, we have extended the powerful frequency-domain picture of femtosecond 
two-photon absorption to a regime of considerable absorption yields, exceeding the weak- 
field yields by more than two orders of magnitude. It corresponds to intermediate field 
strengths, where the interfering pathways are the weak-field (non-resonant) pathways of two 
absorbed photons as well as additional four-photon pathways of three absorbed photons and 
one emitted photon. The picture is based on 4*^-order perturbation theory, which includes 
both the 2"'^ and 4*^* orders associated, respectively, with the two- and four-photon pathways. 
The relative contribution of the 4*''-order absorption amplitude increases as the field strength 
(pulse intensity) increases. This intermediate-field regime is distinguished from the strong- 
field regime where no perturbative description is valid. The strong-field regime is actually 
the one that all the other past multiphoton control studies, which have deviated from the 
weak-field regime, have focused on [TTl UHl [191 ISO]- Our previous work [16] has focused on 
the family of spectral phase patterns that are anti-symmetric around half the two-photon 
transition frequency {ujfg/2). We have found this family to enhance the intermediate-field 
two-photon absorption relative to the unshaped transform-limited pulse, when the central 
spectral frequency is properly detuned, to the red or to the blue (depending on the system), 
from ujfg/2. 

In the present work we systematically study in detail the intermediate-field coherent con- 
trol and the corresponding interference mechanisms, including their dependence on the pulse 
spectrum and its detuning from ujfg/2. The absorption dynamics in the intermediate-field 
regime is compared with the one in the weak-field regime. The model system is the sodium 
(Na) atom. As a test case for femtosecond phase control the study uses the family of shaped 
pulses having a vr spectral phase step, which in the weak-field regime allows high degree of 
control over the full accessible range of the non-resonant two-photon absorption. Section [IT] 
presents and elaborates on the extended frequency-domain 4*'*-order perturbative theoret- 



ical description. The Na intermediate-field control results are presented in Sec. |III[ They 
include experimental results, exact non-perturbative results calculated by the numerical 
propagation of the time-dependent Schrodinger equation, and perturbative results calcu- 
lated numerically using the frequency-domain 4*'*-order formulation. The formers are used 
to validate the latter. Then, in Sec. [IV| the perturbative results are analyzed and discussed 
based on their corresponding frequency-domain description, which allows the identification 
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of the interference mechanisms leading to the different intermediate-field features. 



II. INTERMEDIATE-FIELD THEORETICAL DESCRIPTION 

The atomic femtosecond two-photon absorption process we consider is from an initial 
ground state l^f) to a final excited state |/), which are coupled via a manifold states \n) 
having the proper symmetry. The spectrum of the pulse is such that all the \g)-\n) and 
couplings are non-resonant, i.e., the spectral amplitude at all the corresponding 
transition frequencies is zero: \E{ujgn)\ = \E{u!fn)\ = 0, except for the |/)-|^r) resonant 
coupling for which in general \E{u!fnr)\ 7^ 0. The corresponding excitation scheme is shown 
schematically in Fig. [Tj 

Within the present intermediate-field regime, the time-dependent (complex) amplitude 
af{t) of state |/) at time t, following irradiation with a (shaped) temporal electric field B{t), 
can be validly described by 4*^-order time-dependent perturbation theory. So, in general, it 
includes non-vanishing contributions from both the 2"*^ and 4*^^ perturbative orders: 

afit) = af\t) + af\t), (1) 

with 

af\t) = -—^Hfm/J-^g / e{ti)e{t2)exp[i{ujf^ti + uj^gt2)]dtidt2 , (2) 

1^ „ J -00 J~oo 

af\t) = - fJ'fkfJ'klfJ'lmfJ'mg / / / / ^('^l)^^(^2)^^(t3)^^(i4) 
k,l,m J ~oo J -00 J -00 J ~oo 

X exp [-i{ujfkti + LUkit2 + iOlmh + i^mgU)] dtidt2dt:idti , (3) 

where = {i\fi\j) is the dipole matrix element between a pair of states and Uij = 
{Ei — Ej)/h is the corresponding transition frequency. The 2"'^-order term aj\t) by itself cor- 
responds to the weak-field regime. The intermediate-field final amplitude Af = af(t ^ 00) 
of state I/) after the pulse is over (i.e., t 00) can be expressed as 

A, = A!f + ^K (4) 

with A^j"* = a^p{t 00) and A^^^ = a^^\t 00). The final |/) population Pf = \Aff = 
[9ft(A/)]^ + [$5(y4/)]^ of state |/) reflects the degree of two-photon absorption. 

This perturbative description allows a transformation into a frequency-domain pic- 
ture. Within the frequency-domain framework, the spectral field of the pulse E{oo) = 



|£'(a;)| exp is given as the Fourier transform of B{t), with |-E'(co')| and ^{u) being, 

respectively, the spectral amplitude and phase of frequency u. For the unshaped transform- 
limited (TL) pulse, ^{uj) = for any u. We also define the normalized spectral field 

E{uj) = E{uj)/ \Eo\ 



Eiuo] 



exp [i$(co')] that represents the pulse shape, where \Eq\ is the 
peak spectral amplitude. This allows to clearly distinguish in the expressions given below 
between the dependence on the pulse intensity and the dependence on the pulse shape. The 
maximal spectral intensity Jo is proportional to |-EoP (-^o oc I-EqP). Different values of Jq 
correspond to different temporal peak intensities Jtl of the transform-limited (TL) pulse. 
As shown before for the weak- field regime [B], the 2'^'^-order amplitude A^J"* is given by 



4(2) 



1 

ih? 



Eof A^'Kujf,) , 

DO 

E{uj)E{n - uj)duj 



(5) 
(6) 



where Ufg is the \g)-\f) transition frequency and fijg is the corresponding real effective non- 
resonant two-photon coupling. It is given by fi'j^ = ^'''"-^o ' "^i^h c^o being the carrier 
frequency of the pulse. Eqs. (5)-(6) reflect the fact that Aj- coherently interferes all the 
non-resonant two-photon pathways from l^f) to |/) of any combination of two absorbed 
photons with frequencies u and u' = Ufg — u. Several such two-photon pathways are 
shown schematically in Fig. [T] The phase associated with each two-photon pathway is 
^{uj) + ^{ufg — u). So, with the TL pulse all these pathways acquire zero phase and thus 
interfere one with the other in a fully constructive way. With a given spectrum |£'(ci;)|, this 



leads to the maximal 



.(2) 



and the maximal weak-field non-resonant two-photon absorption. 
The 4*'^-order amplitude term A^^^ is much more complicated than A^^ and we have 



calculated it to be given by 



A 



(4) 



1 



En 



1 (4)on-rcs < (4)ncar-res 



^(4)on-res ^ ,^^(2) (Q) 



A 



(4)near-res 
/ 



d5-^A^^\ujfg-5)A^''\5) 



(7) 
(8) 
(9) 
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where A^'^\n) is defined in Eq. and 



/oo 
E{uj + An)E*{uj)diu 
-oo 



d5'-E{uJfn. +An- 6')E*{Ufn^ - 5') 



(10) 

(11) 

(12) 



This set of equations reflects the fact that Aj interferes all the four-photon pathways from 
\g) to I/) of any combination of three absorbed photons and one emitted photon. Several 
typical four-photon pathways are shown schematically in Fig. [T] 

Each four-photon pathway can actually be divided into two two-photon parts: (i) a 
non-resonant transition of two absorbed photons u and u' with a frequency sum oi Q = 
uj + u' = ujfg — 6, and (ii) a Raman transition of two photons ur and u'p^ with a frequency 
difference of AQ = uor — u'j^. The border line between these parts is detuned by 6 from 
either |/) or \g) according to whether, respectively, part (i) or part (ii) comes first (see 
Fig. 1). The and ^Wnear-res |-gj,j^g ^(4) ^j^qg^ (7)_(9)j interfere, respectively, these 

on-resonant {6 = 0) and near-resonant {6 ^ 0) four-photon pathways. The on-resonant 
pathways are excluded from Iqj ^]-^g Cauchy's principle value operator p. The 

integration taking place in these terms over the corresponding pathways is expressed using 
the product of two parameterized amplitudes, A^'^\^l) and A^^\A^l), which originate from 
the different two-photon parts of the four-photon pathways. The amplitude A^'^\n) interferes 
all the non-resonant two-photon transitions [parts (i) above] with transition frequency fl, 
while the amplitude A^^\AQ,) interferes all the Raman transitions [parts (ii) above] with 



transition frequency A^l. As given by Eqs. (10)-(12), A^^\AQ) includes two components. 
The first component is y^C^ion-res^?) interfering all the Raman transitions that are of non- 
resonant nature, with /i^^ and fi'jj being the \g)-\g) and |/)-|/) real effective non-resonant 
Raman couplings due to all the non-resonantly coupled states \n). They are given by fi 
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. Hence, their sum 



{Hgg + jiff) [appearing in Eq. (11)] is either positive or negative, depending on the physical 



system and pulse spectrum. The second component of A'^^\AVt) is A^'^^^^'^ interfering all the 
Raman transitions that are of resonance-mediated nature via In^.), with 5' being the detuning 
from \nr) (see Fig. [T]). These on-resonant [5' = 0) and near-resonant {5' ^ 0) transitions are 



interfered separately, respectively, in the first and second terms of A^^'^^^'^ [Eq. (12)]. The 
Cauchy's principle value operator p excludes the on-resonant transitions from the second 
term. 

For a given physical system and a given pulse shape E{uj), a non-zero Aj is proportional 
to I-EqP while a non-zero A^^^ is proportional to I-EqI^, i-e., their ratio is proportional to l-^oP 
or, equivalently, to Jq. For a given Jq, the relative magnitude and relative sign between the 
real parts 3?[y4j^''] and 3fJ[Aj^''] are generally determined by the pulse shape E{uj) and by the 
magnitudes and signs of the different Raman couplings (/igg, fij-f, and l/i/nrP)- The same 
applies also for the imaginary parts 53 [A ^^^] and ^^[A^'^'*]. 

In the present work, for a set of intensities Jq, the final |/) population Pj (i.e., the degree 
of two photon absorption) is controlled via the pulse shape E{uj). Three different spectra 
E{uj) are being considered, with the control knobs being the various spectral phases 
The three spectra are chosen such that their central spectral frequency uoq is of no detuning, 
blue detuning, or red detuning from ujfg/2, with E{u)fg/2) ^ 0.5 for the detuned cases. 
This set of spectra generally corresponds to the typical case, where A'^^^ is negligible relative 
to A^'^'> in the weak-field limit and becomes comparable to A*^^^ in the upper intermediate-field 
limit. These spectral shifts do not change the sign of the various Raman couplings and hardly 
change their magnitude. They do however affect the relative amplitude associated with the 
different interfering pathways and, thus, their overall interference result. For example, as 
shown below, in the cases studied here the spectral change from red- to blue-detuning leads 
to a change in the relative sign between 53 [A y^''] and ^'[Ay^''], and thus to a change in the 
nature of their interference from destructive to constructive. 

In general, the dynamics and interference mechanisms discussed in this work involve intra- 
term as well as (intensity-dependent) inter-term interferences involving the multiphoton 
pathways of the 2"^^- and 4*'*-order perturbative terms. 



III. RESULTS 



The physical model system of the study is the sodium (Na) atom [23]. It includes the 
3s ground state as \g), the 4s state as |/), the manifold of p states as the \n) manifold, 
and the 7p state as 1^^). So, A^p = A^^} and A^'^^ = A^^] . The transition frequency 
ujfg = u^s^ss = 25740 cm~^ corresponds to two 777-nm photons and the transition frequency 
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ujfn^ = ujjp^4s = 12801 cm~^ corresponds to one 781.2-nm photon. The atomic sodium 
is irradiated with phase-shaped hnearly-polarized femtosecond pulses having an intensity 
spectrum of ~5-nm FWHM-bandwidth (~180 fs TL duration) centered around a wavelength 
Ao tunable between 773 to 780 nm. As a test case, for three different values of the central 
wavelength, the present study uses the family of shaped pulses having vr-step spectral phase 
patterns. In the weak-file regime [6] this family allows high degree of control over the full 
accessible range of the two-photon absorption, from zero to the maximal level (induced 
by the TL pulse). Each such pattern is characterized by the vr-step position Ugtep, with 
$(cj) < Ustep) = -7r/2 and > u^tep) = 7r/2. 

The understanding and analysis of the intermediate-field control mechanism is conducted 
below using the frequency-domain picture given in Eqs. (|4])-( 12 ) and corresponding numerical 
results for the Na system. However, prior to the analysis, the extended perturbative picture 
and results are validated by a comparison to exact non-perturbative results that their own 
validity is confirmed first by a comparison to experiment. The exact non-perturbative re- 
sults have been calculated by the numerical propagation of the time-dependent Schrodinger 
equation equation (TDSE) using the fourth-order Runge-Kutta method. The theoretically 
considered manifold of p-states is from 3p to 8p, including all the (1/2 and 3/2) fine-structure 
states [231. 



A. Non-perturbative calculations vs. experiment 

Experimentally, atomic sodium vapor is produced in a static chamber at 300°C (Na partial 
pressure of ~0.1 Torr) with 10-Torr Ar buffer gas. It is irradiated at a 1-kHz repetition 
rate with shaped femtosecond laser pulses of three different central spectral wavelengths: 
Ao=773, 777, and 779.5 nm. The corresponding spectral intensity bandwidth (FWHM) is, 
respectively, 5.5, 4.5, and 5 nm. The 773-nm and 777-nm spectra are Gaussians, while 
the 779.5-nm spectrum is a modified Gaussian having a slight asymmetry toward short 
wavelengths. The slight change in the experimental spectral shape between the three cases 
results from technical limitations and is of no significance to the analysis and discussion 
presented below. The laser pulses undergo shaping in a 4/ optical setup incorporating a 
pixelated liquid-crystal spatial light phase modulator [21]. The effective spectral shaping 
resolution is Sushaping = 2.05 cm~^ (0.125 nm) per pixel. The experiment is conducted with 
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different pulse energies. Upon focusing, the corresponding temporal peak intensity of the 
transform-limited (TL) pulse at the peak of the spatial beam profile ranges from 

5x10^ to 7x10^° W/cm^. Following the interaction with a pulse, the Na population excited 
to the 4s state radiatively decays to the lower 3p state. The fluorescence emitted in the 
decay of the 3p state to the 3s ground state serves as the relative measure for the excited 
4s population Pf = P4S. It is optically measured at 90° to the beam propagation direction 
using a spectrometer coupled to a time-gated camera system. The measured signal results 
from an integration over the spatial beam profile. 

Figure [2] compares the non-perturbative theoretical results (solid lines) with the exper- 
imental results (squares) for the two-photon absorption in Na. Each column in the figure 
corresponds to a different spectral case of Ao=773, 777, or 779.5 nm with different pulse 
energies, i.e., different ^jj-^g pulse energy increases from top to bottom within 

a single column. The traces show the final 4s population as a function of the vr phase 
step position Ustep- Each of the traces is normalized by the final 4s population P4s,tl excited 
by the corresponding TL pulse. The weak-field vr-traces are given in the first-row panels 
[(a)-panels] of Fig. g The non-perturbative theoretical results shown in Fig. [2] account for 
the experimental integration over the spatial beam profile. Each of the presented traces 
results from an appropriately-weighted integration over a set of calculations conducted each 
with a different single value of Iq. 

As can be seen, there is an excellent agreement between the experimental results ("real 
experiment") and the non-perturbative results ("computer experiment"). Hence, the accu- 
racy of the latter is confirmed for the present intermediate-field excitation of Na. 

B. Perturbative calculations vs. non-perturbative calculations 

Next, the confirmed non-perturbative calculations are used to validate the intermediate- 
field perturbative results and to identify the intensity limit of the intermediate-field regime 
for the present Na excitation, i.e., the intensity up to which no perturbative order beyond 
the 4*'' one is needed to be included. From this point on, the analysis and discussion are 
conducted based on the theoretical results. So, for completeness, the three intensity spectra 
considered theoretically are of a perfect Gaussian shape with 5-nm bandwidth (FWHM) 
centered around Ao=773.5, 777, and 780 nm. As mentioned above, these Aq values correspond 
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to blue detuning, no detuning, and red detuning of ujq from uj4s,3s/'^, with E{u!4s,3s) — 0.5 
for the detuned cases. 

Figure [S] compares the theoretical non-perturbative results (thick gray lines) and pertur- 
bative results calculated numerically using Eqs. (|4])-( 12 ) (thick black lines). The perturbative 
calculations use a frequency grid with a bin size equal to the experimental spectral shaping 
resolution Sushaping- The presented vr-traces are given on a TL-normalized scale of P4s/P4s,tl 
(left-hand y-axis) as well as on an absolute scale of P^s (right-hand y-axis) . The x-axis is the 
71 phase step position Ustep- Also here, each column corresponds to a different spectral case 
of Aq with different (single-valued) intensities Iq, i.e., different /tl- The Iq increases from top 
to bottom within a single column. Shown are several examples out of the full set of results. 
The weak-field traces are given in the first-row panels [(a)-panels] of Fig. [s] Actually, on a 
TL-normalized scale they are all identical and independent of Aq [6J. For comparison, the 
TL-normalized weak-field trace is also given in all the other panels (thin black line). As 
can be seen from the figure, the perturbative results reproduce the exact non-perturbative 
results up to the Jq that corresponds to a TL peak intensity of /tl = 2.5 x 10^° W/cm^ 
[(c)-panels of Fig. [sj. This is the intensity limit of the 4*^-order intermediate-field regime for 
the present Na excitation. As the last-row panels [(d)-panels] of Fig. [s] show, the 4*'^-order 
perturbative description is not sufficient above this intensity. The corresponding intensity 
limit of the weak-field regime, where Af ^ , is /tl ~ 5 x 10® W/cm^. 



C. Intermediate-field features of the vr-trace 



The prominent features of the intermediate-field vr-traces are presented below. As pre- 
viously shown |6] and can be seen in Fig. [3] [the (a)-panels], the weak- field TL-normalized 
TT-trace is symmetric around ujstep = '^4s,3s/2 (777 nm) and is identical for any uq (Aq). 
Its shape is determined only by the spectral bandwidth of the pulse. Also, since it is TL- 
normalized, it is independent of the intensity Iq [see Eq. ([s])]. However, when deviating from 
the weak-field regime with the 4*^* perturbative order playing a role, the TL-normalized vr- 
trace losses its weak- field symmetry and its shape becomes dependent on both ujq and Iq. 
At a given uq, the degree of deviation from the weak-field shape depends on Iq. 

As mentioned above, in the weak-field regime the maximal non- resonant two-photon 
absorption is induced by the TL pulse. Additionally, as can be seen in the figure, the same 
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maximal weak-field two-photon absorption is also induced by the shaped pulse with oOstep = 

l^4s,3s/2 (777 nm), i.e. -P4s,step@aj4s,3s/2 = -Rls.TL or -Rstep@a;4s,3s/2 = P4s,step@oJ4s,3s/2/ P4:S,TL = 1- 

Both pulses induce fully constructive interferences (i.e., zero relative phase) between all 
the \g)-\f) two-photon pathways. Conversely, in the intermediate-field regime a vr-phase 
step at Ustep = "^43,35/2 leads to a two-photon absorption that, in general, is different from 
the TL absorption. The difference is highly pronounced when uq is detuned from uJis,3s/'^ 
[Fig. [3]- columns(l) and (3)] and is weakly pronounced when coq = (^4^ 3^/2 [Fig. [3]- column 
(2)], with -R4s,step@a;4s,3s/2 < 1 for a bluc detuning, -R4s,step@a;4s,3s/2 ~ 1 for no detuning, 
and -R4s,stcp@a;4s 3^/2 > 1 for a red detuning. For example, at the intermediate-field limit, 
with Jo of Jtl = 2.5 x 10^° W/cm^, Rstep®uj4,s,3s/'2 reaches a value of 0.65, 1.15, 2.0 for 
Ao=773.5 nm [Fig.glc)], Ao=777 nm [Fig. |3|(2c)], and Ao=780 nm [Fig.g3c)], respectively. 
The corresponding degree of attenuation or enhancement over the TL absorption , i.e., 
|i?step@a;4a,3s/2| , iucrcascs as Jo increases. Also, as Jq increases, this enhancement /attenuation 
effect occurs (with variable corresponding values) over an increased range of ujgtep around 

'^4s,3s/2. 

In the spectrally detuned cases of ^ uji^sfisl'^ with £'(0745,3^/2) ~ 0.5 (as is the case 
here; see above) a vr phase step that is positioned around c^o, i-e., ujstep ~ i^o, leads to 
intermediate-field values of the TL-normalized trace [Fig. [3]- columns (1) and (3), thick 
solid lines] that are systematically either attenuated or enhanced relative to the weak-field 
trace [Fig.|3]- columns (1) and (3), thin solid lines]. The attenuation occurs for blue-detuned 
Ao and the enhancement occurs for red-detuned Aq. Also here, the degree of attenuation 
or enhancement relative to the weak-field case at a given io^tep increases as Jq increases. In 
the red-detuned case, at high enough intensity Jq, the corresponding two-photon absorption 
even exceeds the TL absorption [see Fig. [3|^3c)]. 

An additional intermediate-field feature is the dip observed when uJstep = '^ip^ 
(781.2 nm), which does not exist in the weak-field trace. It occurs when the field amplitude 
E{ujfn,.) = E{uj7p^4s) at LJip^4^s significantly deviates from zero. Hence, this 781.2-nm dip 
is very prominent for Ao=780 nm [Fig. |3]- column (3)], weakly noticeable for Ao=777 nm 
[Fig. [3]- column (2)], and hardly noticeable for Ao=773.5 nm [Fig. [3]- column (1)]. 

Lastly, it is important to mention that within the intermediate-field regime, due to the 
4*^-order term ^4'^'', the absolute 4s population P4s,tl excited by a given shaped pulse 
strongly deviates from the weak-field intensity dependence of (Jq)^ in all the Aq cases. Specif- 
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ically, it applies also to the non-detuned case of Ao=777 nm, even though the corresponding 
intermediate-field shape of the TL-normalized vr-trace does not deviate much from the weak- 
field shape. This can be seen in Fig. |3| for example, in the absolute population values that 
are induced by the TL pulse (i.e., the asymptotes of the traces) at different intensities. 



IV. DISCUSSION 



A. General considerations 



The discussion below analyzes the coherent source for the intermediate-field features 
described above, using the frequency- domain picture given in Eqs. (|4|)-(|12|). The different 
intermediate- field features, i.e., the deviation of the intermediate- field trace from the weak- 
field trace, originate from the different dependence of the 2"''- and 4*^-order amplitudes, A^'J 
and A^^] , on the pulse shape. 

For a given spectrum the on-resonant term ^4'^''°'^''^*''^ of ^4'^'' [Eq. (js)] is propor- 

tional to A^^] . It is so since A^^\AQ = 0) depends only on the spectral intensity: all 

the corresponding Raman transitions involve two identical photons and their amplitudes in 

J I J I _ _ _ 2 

Eqs. (10)-(12) are of the form E{uj)E*{uj) = E{uj) . Thus, the difference in the phase 



dependence of A^^] and ^4^^ originates only from the near- resonant term ^W^ear-res ^4"^'*, 
which interferes all the near-resonant four-photon pathways of non-zero detunings 5 7^ 
from 3s or 4s (see above). 

As Eq. (9) shows, Af is given by a proper p-integration over all the non-zero 
5 values. Its phase dependence originates from the phase dependence of the integrands 

^(2) 

('^4s,3s~'^) and A'^^^S). Their dependence on the spectral phase pattern ^{u) is reflected 
in their values for a given 5 and, more importantly, in their functional dependence on 
5. So, the p-integration yields different results for different phase patterns. Due to the 
1/6 weighting (and its sign change for ±|5|), the p- integration result is dominated by the 
integration over small values of |5| and is highly sensitive to the degree of symmetry of the 
integrand A^'^^uj/^g^'^s — 5)A^^\5) around 5 = 0, i.e., how different are its values for ±|5|. 
Below, the intermediate-field two-photon absorption is analyzed based on the 5-dependence 
of A(2)( 

^As,'is ~ a-iid A^^^{5) for several representative cases of the vr-step position ujgtep 
with the different Aq. 
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The analysis of the amphtude A^^\5), which is contributed by all the Raman parts of 
the 5-detuned four-photon pathways, is simplified by including its component A^'^'^^^^ (6) [see 



Eq. (10)] only for the study of the dip feature at oOstep = i^7p,4s (781.2 nm). As described 
above, A^'''^'^'^) (5) interferes those Raman parts that are resonance-mediated via the 7p state. 

This line of analysis is supported by the discussion given below and by the perturbative 
results presented in Fig. [4] for the vr-traces at Iq of Jtl = 2.5 x 10^° W/cm^ (i.e., the 
intermediate- field limit). In addition to the real TL-normalized vr-traces (gray thick lines), 
which are also shown in the (c)-panels of Fig. [sj the figure displays the vr-traces (black 
thick lines) calculated with A^^\6) = A^''°''-'^''^\6) , i.e., with artificially setting A^'''''^^{S) 
to zero for any 6. As can be seen, the resonance-mediated Raman transitions via 7p are 
of significance only for discussing the dip at ujgtep = "-^7^,45 (781.2 nm). For completeness, 

the figure also displays the weak-field trace (black thin lines; shown also in Fig. [3]) that 

(2) 

originates only from Ar and, thus, is not affected by changes in the Raman part. 



The non- resonant Raman term A^^°^~'^'^^^\5) [Eq. (11)] coherently integrates all the non- 
resonant Raman amplitudes E{u! + 6)E*{u!) contributed by all the possible pairs of photons 
with a frequency difference of 6 {00 is scanned across the spectrum). With ^{uj) = (TL 
pulse) or any vr spectral phase step, the amplitude contributed by any such pair of photons 
is a real (positive or negative) quantity and, thus, so is the resulting A^^°^~'^'^^^\5) . So, with 
these phase patterns, the relation A('^°"-'''="^)(-|5|) = A^''°''-'''^^\+\6\) holds for any 6, i.e., 
^(non-resfi) j^^-j is Symmetric around 5 = 0. This symmetry can be seen in Fig. [sjthat shows, as 
an example, the value of A^^°'^~'^'^^^\5) (gray thick lines) as a function of 5 for the TL pulse 
and for the shaped pulses with Ugtep = i^4s,3s/2 and Ust^p = uq. Each column corresponds 
to a different Aq (cuq)- The x-axis values are actually the normalized detuning values 6/Auj, 
with Auj being the bandwidth of the intensity spectrum. As can also be seen in the figure, 
^(non-rcsfi) j^^-j is maximal at 5 = 0, with a value that depends only on the spectral intensity 
E{uj) , i.e., is independent of $(a;). 

The other term A^'^\lu4s,3s — [Eq. (j6|] coherently integrates all the non-resonant two- 
photon transition amplitudes E{uj)E{lj4s^3s — contributed by all the possible pairs of 
photons with a frequency sum of uj4s,3s — Also here, with $(c<j) = (TL pulse) or 
any vr spectral phase step, the amplitude contributed by any such pair of photons is a real 
(positive or negative) quantity and, thus, so is the resulting A^'^\uj4s,3s — S). Figure [s] displays, 
together with the Raman part data, the corresponding A^'^^ (t<^4s,3s — S) (black thick lines) as 
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a function of S/Au in the different cases. The zone of small \S\ around S = 0, which is the 
most contributing to j\^^'>^^^^~^'^'^ (^ggg above), is also indicated schematically. As can be seen 
from the figure, the 5-dependence of A^'^^cUis^^s ~ 5) around 5 = is different from one case 
to the other and depends on both the spectral phase pattern and Aq {ujq). Each panel also 
shows, for comparison, the A^'^\uj4s,3s — S) trace of the corresponding TL pulse (black thin 
line). 



B. Intermediate-field two-photon absorption: shaped pulse with ujstep = '^As,3s/'^ 
vs. TL pulse 



The first-row panels [(a)-panels] of Fig. |5] correspond to the TL pulse of the different Aq 
(ujq) cases. Since the TL pulse induces fully constructive interferences among all the two- 
photon pathways that contribute to A^'^\Q) for any Q, the corresponding value A^^l^ = 
i^4s,3s — is actually the maximal one (positive and real) for any given 6. As can be seen 
from Eq. (6j), A}^l{Q = uj^s^-^s — is actually a convolution of the corresponding spectrum 
£■(0;) . Thus, with a Gaussian spectrum around ujq, it is peaked at fipeak = '^^0, i-e., at 
Speak = i^4s,3s — Scuq- So, for the different cases of Uq one obtains the following behavior: (i) 
For luq > LU4s,3s/'^ [Ao =773.5 nm. Fig. [sjl^la)] - 5peak < and A^'^\uj4s,3s — monotonically 
decreases around 6 = as 6 increases from negative to positive values, i.e., A^'^\uj4s,3s — \S\) > 
A^'^\uj4s,3s + for small \6\; (ii) For uo < uj4s,3s/2 [Aq =780 nm. Fig. [s^Sa)] - Vak > 
and A^'^\uj4s,3s — 5) monotonically increases around 5 = upon the negative-to-positive 
increase of S, i.e., A*^^) (c<j4s 3^ — \5\) < A^'^\uj4s,3s + l*^!) for small \5\; (iii) For cuq = uj^sfis/'^ 
[Ao =777 nm. Fig. [5|^2a)] - 5peak = and A^'^\ijj4^s,3s — 5) is symmetric around (5 = 0, i.e., 

A^'\uO,s,3s-m = A^^K^As,3s + m- 

The second-row panels [(b)-panels] of Fig. |5] correspond to the shaped pulse with uJstep = 
i^4s,3s/2 (777 nm). Generally, when a vr phase step is positioned at Ustep, the value of 
^(2)(^ = 2ujstep) is equal to the corresponding TL value Axl(^ — '^^step) due to fully 
constructive interferences among all the involving two-photon pathways. As for the TL 
pulse, the phases associated with these pathways are all zero. However, as deviates from 
2ujstep^ the value of A^'^\VL) gradually reduces with comparable magnitude for positive and 
negative deviations. Hence, a vr-step at uJstep = i^4s,3s/2 yields a peak of A^'^\uj4s,3s — 5) at 
6 = (with the TL value). When ujq = a;4s,3s/2 [Ao=777 nm; Fig. |§2b)], this peak IS a 
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global one and A^'^^to^s^^s — ^) has a perfect symmetry around 5 = 0, i.e., A^'^\ijjis^'is — = 
A^^\iU4,^3s + When coq ^ cu4s,3s/2 [Aq =773.5 and 780 nm; Fig. [s^lb) and (3b)], this 
peak is a local one with only an approximate symmetry of A^'^\u!4^s,3s — S) around 5 = 0, i.e., 
A^^\uj^s,3s - ^ A(2)(^4,^3, + \6\) for small \6\. 

Based on the above analysis of A^^°^~'^^^^\5) and A^'^\uJis,3s — 5), considering also the 
integrand factor 1/6 that is anti-symmetric around 5 = 0, the magnitude of ^Wnear-res 



can be compared between the TL pulse case 



A 



(4)near-res 
4s,TL 



and the shaped pulse case of 



^step — l^4s,3s/2 

A (4) near-res _ 
^4s,step@a)4s,3s/2 



^(4) near-res 

4s,step@t<J43^3s/2 



i) In the no-detuning case of Ao=777 nm - A^^J^^^'^*^^ = 
0, since the amplitudes contributed to the p-integral by the four-photon 
pathways of positive and negative detuning ±|5| are equal and, thus, cancel out each 
other; (ii) In the red- and blue-detuning cases of Ao=773.5 and 780 nm 



A 



(4)near-res 
4s,TL 



> 



A 



(4)near-res 



since the amplitudes contributed by the four-photon pathways of detun- 
ings ±|5| are approximately equal one to the other for the shaped pulse of ujstep = i^4s,3s/2 
while they are significantly different one from the other for the TL pulse. 

As discussed below, important is also the sign of ^Wncar-res j,gjg^^jyg A^^J , which is 
determined by the pulse spectrum and by the non-resonant Raman couplings sum {fil^ ^s + 
As)- The sign of A^^]^^"^ ^'^'^ relative to ^4^''tl ^he sign of A 



(4) near-res 
4s,step@aJ4s^3s/2 



relative to 



A 



(2) 

4s,step@t<J4s ^33 /2 



are obtained here to be the same for any Aq. However, upon a blue-to-red 



spectral shift of Aq (i.e., from 773.5 to 780 nm) these signs change from positive to negative, 
with the sign flip occurring when passing via Ao=777 nm. The sign change results from a 
change in the relative magnitude of the amplitudes contributed by the negatively-detuned 



and positively-detuned four-photon pathways, while {fi 



3s, 3s 



/i|^4g) keeps its sign and is 



effectively constant over the whole Aq range considered here. 

For comparing the two-photon absorption induced by the TL pulse vs. the absorption 
induced by the shaped pulse of Ustep = i^4s,3s/2, one needs to consider the coherent amplitude 
addition of ^4^-' and A^^] . For the TL pulse and for any shaped pulse with a vr-step phase 



(2) 

pattern, A\J [Eq. (5)] is an imaginary quantity, i.e., A 
quantity. When only ^{°on-res_R) -g included in A'^^^ (see above), the corresponding ^Wnear-res 



(2) 
4s 



53 



A 



(2) 
4s 



while 744^-* is a complex 



component is real while the corresponding A 

1(4) 



(4)on-res 
4s 



component is imaginary. Due to the 1/i 



preceding factor of ^4^ [Eq. (7)], they contribute, respectively, solely to the imaginary and 



real parts of A_ 



(4) 
4s ) 



i.e., 53 



A 



(4) 
4s 



oc4y"™and3fJ 



A 



(4) 
4s 



oc A^^g°^~'^'^^ . Thus, the interferences 



between the pathways groups of A\J and of AlJ actually take place between those included 
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in A^^} and those included in 53 



A 



A (4)near-res 

I.e., m A\J 



4s 

Since the following relations are satisfied here: (i) A\ . ^ ,„ 



I (2) 



a{2) ^ A2) .... 
^4s,TL — ^4s,*5 lllj 



Pis 
I (2) 



1^ 



4s 



A 



(2) 



4s 
(4) 



53 



A 



3? 



A 



(4) 
4s 



fiii) 3? 



A 



(4) 
4s 



(2) /Ik 

oc ^^43 [see Eq. (8)], and (iv 



^44^"^ oc Jo while ^4^'^ oc Jq [see Eqs. (5) and (7)], the TL-normalized vr-trace value -Rstep@w4s,3s/2 
(real and positive) corresponding to uj^tep = ^^4s,3s/2 is given as 



R 



'Step&uJ4s,3s/2 



4s ,step@[i;4a _ 3s /2 



4s,TL 



(13) 



l + 5> 


.(4) 

4s,step@ci;4i,_3s/2 


mIsL 


2 

+ 




.(4) 
4s, step 




l + $5 


.(4) 
^4s,TL 


mIsL 


2 

+ 


5? 


.(4) 
^4s,TL 



(2) 
4s,* 



mIs!* 



I I Z I I 

|1 + -^step@LJ4s,3s/2-^o| + I'^-^ol 



|l + i^TL/o| +|/t/o| 

where the pulse shape dependence enters only via the K factors. Accounting for the mag- 



nitude of 



A 



(4)near-res 
4s,TL 



VS. 



A 



and for the sign of relative to ^44^'*, 

one obtains the following behavior for the different Aq cases: (i) For blue detuning of 



(4)ncar-rcs 
4s,step@aJ4s,33/2 



^0 > ^4s,3s/2 (Ao=773.5 nm): K 



step(aLiJ4s_3 



,/2 > 0, i^-TL > 0, \K 



step@ciJ4s,3s/2 I 



TL 



thus R. 



■step@(^4s,3s/2 



< 1; (ii) For no detuning of Uq = c<j4s 3s/2 (Ao=777 nm): 



step@a;4s,3s/2 



and 
0, 



i^'xL = 0, and thus -Rstep@a;4s,3s/2 = 1; (iii) For red detuning of < ll!4^s,3s/'^ (Ao=780 nm): 

Kstep@uJ4s,3s/2 < 0, Kn[L < 0, | -ft'step@tJ4s,3s/2 1 < I-^^Tl], and thuS -Rstep@a;4s,3s/2 > 1- In 

other words, the two-photon absorption induced by the shaped pulse of Ustep = i^4s,3s/2 
is lower, equal, or higher than the TL absorption according to whether ujq (Aq) is blue- 
detuned, non-detuned, or red-detuned from 0745^3^/2, respectively. One can also see that, 
when |-Rstep@a;4s 3^/2! 7^ I5 the value of |-Rstep@a;4s 3^/2! increases as Jq increases. This entire 
intermediate-field behavior is indeed the one observed in the results of Fig. [3} except for a 
small deviation of Rstep@uj4s,3s/2 from a value of one in the case of Ao=777 nm [for example, in 
Fig. |3|^2c) it reaches a value of 1.15]. This deviation originates from the resonance-mediated 
Raman term A^^*^^^^ that is excluded from A^^^ in this part of the analysis and is not sym- 
metric around 5 = [as A^'^"""''^'^^) (5) is]. 

C. Intermediate-field two-photon absorption: shaped pulse with uJstep ~ 1^0 



The third-row panels [panels (Ic) and (3c)] of Fig. |5] correspond to the shaped pulse 
with LOgtep = "^o in the detuned Aq cases of 773.5 nm and 780 nm (the corresponding case 
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with Ao=777 nm is actually the one already considered in Sec. IV B). It is considered here 



as the representative case for the region of ujstep ~ i^o; where a systematic attenuation or 
enhancement of the intermediate-field TL-normalized absorption relative to the weak-field 
TL-normalized absorption occur [Figs. [3|lc) and (3c), Figs. |4]^a) and (c)]. The attenua- 
tion or enhancement correspond, respectively, to the blue (Ao=773.5 nm) or red detuning 
(Ao=780 nm) of ujq from (UJ4s,3s/2. Detailed analysis of the present conducted above 

for the TL pulse and for the shaped pulse with Ugtep = "^45,35/2, leads to the conclusion that 

the only difference here is that 4!lstep@<^o ^ ^Stl while ^i5step@^4.,3./2 = ^2tl- All the 
other qualitative conclusions regarding the different amplitude components in terms of their 
magnitude, sign, and real/imaginary character are exactly the same for both the Ustep = ^0 
and Ustep = t^4s,3s/2 pulses. 

In general, the ratio RRstep between the intermediate-field and weak-field TL-normalized 
absorption corresponding to the shaped pulse of a given Ugtep is given by 



RR. 



'Step@aJstep 



pinterm-field 
pweak-field 

\tep@UJstep 



(14) 



1 + Q 



4s,step@ 



(2) 

4s,Step@tJstep 



+ 



3? 



A 



(4) 

4s,step@aJst, 



(2) 

4s,Step@UJstep 



1 + ^ 



A 



(4) 

4s,TL 



(2) 

4s,TL 



+ 



.(4) 
^4s,TL 



(2) 

4s,TL 



|l + Kstep@oJstep^o\ + I'^-^ol 



where the pulse shape dependence enters only via the K factors. As can be seen, there 



is a close similarity between this equation and Eq. (13). So, based on the above analysis 
conclusions, one gets the observed intermediate-field behavior for Ugtep = ^o- (i) For blue 
detuning of t^o > uj4s,3s/'^ (Ao=773.5 nm): K, 



stcpOo^o 



> 0, i^TL > 0, \K 



stepOo^o I 



TL 



and thus -R-RstepOojo < 1 (attenuation); (ii) For red detuning of Uq < U4s,3s/'^ (Ao=780 nm): 



K 



step@ui() 



< 0, i^TL < 0, \K 



stepOojo 



< I-R'tlI, and thus RR. 



'Step@a;o 



> 1 (enhancement). Also, 
for both Ao values, the value of |-R-Rstep@a;ol increases as Iq increases. In the red detuning 
case (Ao=780 nm), this increase leads eventually to an intermediate-field absorption that 
exceeds the intermediate-field TL absorption [Fig. [3]^3c)]. It worth emphasizing that the 
intermediate-field attenuation/enahn cement effect considered here is relative to the weak- 
field case of the same shaped pulse, while for the pulse with Ugtep = i^4s,3s/2 it is considered 
above relative to the intermediate-field TL absorption. 
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D. Intermediate-field two-photon absorption: shaped pulse with 00step='^7p,As 



As discussed above with regard to Fig. |4| the intermediate-field dip feature at Ugtep = 
^7p,4:s (781.2 nm) originates from the inclusion of the Raman term A^^^^^"* in the complete 
Raman term A^^'^ [Eq. (10)-(12)]. As described above, A'^^^^^'> interferes only those Raman 
parts that are resonance-mediated via Ip. Each of them is either on resonance or near 
resonance with 7p {5' is the corresponding detuning). So, overall, each corresponding four- 
photon pathway is either on- or near-resonance with 4s and either on- or near-resonance 



with 7p. The higher is the field amplitude 



E{ujjp^As) at ^7p,4s the more prominent is the 



dip. So, here, it is most prominent for Ao=780 nm [Fig. |4]^c)]. 

With the inclusion of A^^'^^^\ the terms a^^}^^'^'^'^^^ and A^^}^^''^'^^ become complex quan- 
tities, and thus they both contribute to both the imaginary and real parts of A^^} . So, the 
analysis of the corresponding interference mechanism is much more complicated as compared 
to the above, when A^'^^^^'^ is excluded. Essentially, the interference mechanism leading to 
the dip is very similar to the one we recently identified in weak-field resonance-mediated 
(2-1-1) three-photon absorption [14J, where a vr-step at ^77^ 4^ leads to a strong enhancement 
in the population transfer to the 7p state. Here, for example, a vr-step at uj^p^^s leads to a 
constructive add-up within the p-integral of A^^]'^'^^'^'^'^^ [Eq. (|9)] between the amplitudes con- 
tributed by the four-photon pathways on-resonant with 7p that are of positive detuning + \5\ 
and of negative detuning —\5\ from 4s. The corresponding Raman transition involves the 
absorption of the photon ujjp^^s^^ and the emission of the photon Uj^^^^s. Since E{ujjp^A^s + \5\) 
and E{u7p^4s — are of opposite signs for the shaped pulse of Ustep = ^7p,4s, the sign of the 
detuning 6 (and of the integrand factor 1/6) becomes correlated with the sign of E{ujp^4^s+6), 
leading to the constructive add-up and resulting amplitude enhancement. For the TL pulse, 
the add-up is destructive since the signs of E{ujp^4^s — and E{uip4s + \6\) are the same. 
Upon a detailed analysis, accounting also for the 5-dependence of A^'^^tOis^^g ~ \^\) when 
uJstep = i^7p,4s (781.2 nm) (very similar to the one when \step=780 nm [Fig. |5](3c)]), one 
obtains that, for the present Na excitation, the sign of the amplitude contributed to ^4^^ by 
the four-photon pathways with non-resonant Raman parts [i.e., corresponding to A^^'^^^^] is 
opposite to the sign of the amplitude contributed by those with resonance-mediated Raman 
parts [i.e., corresponding to A^"^^^-^^]. This leads to a dip feature rather than a peak feature. 
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V. CONCLUSIONS 



In conclusion, coherent control of femtosecond two-photon absorption in the intermediate- 
field regime in analyzed in detail using a powerful frequency domain description that is 
based on 4*''-order perturbation theory. The two-photon absorption is coherently induced 
by non-resonant two-photon transitions as well as by four-photon transitions that introduce 
a resonance-mediated nature to the excitation. Their relative contributions to the total 
absorption amplitude depend on the field strength. The corresponding interference mech- 
anism is identified to include intra-group and inter-group interferences involving these two 
groups of multiphoton transitions. The inter-group interferences lead a difference between 
the intermediate-field and weak-field absorption dynamics. Their constructive/destructive 
nature is found to depend on the detuning direction of the pulse spectrum from half the two- 
photon transition frequency; It changes upon a red-to-blue detuning change. The extended 
frequency-domain description and its detailed understanding serve as a basis for femtosecond 
control with rationally-shaped pulses in a regime of significant absorption yields, reaching 
population transfer in the range of 10-40% (depending on the specific system and excitation 
scheme). They also serve as a basis for future extensions to molecular systems, to other 
types of multiphoton processes, and to more complicated excitation schemes. 
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FIG. 1: Excitation scheme of femtosecond two-photon absorption in the intermediate-field regime. 
The indicated levels correspond to the Na atom (not to scale). Shown are pathway examples of 
non-resonant two-photon transitions and four-photon transitions from |/) = 4s to 1^) = 3s. The 
four-photon transitions involve three absorbed photons and one emitted photon in any possible 
order, and thus can be decomposed into two parts: a non-resonant two-photon transition and a 
Raman transition. The border line between these two parts can be either on-resonance or near- 
resonance with 3s or 4s (with detuning 6). The Raman transition itself can be non-resonant due 
to the np states (except for 7p) or on/near-resonance with 7p (with detuning S'). 
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FIG. 2: Experimental results (squares) and non-perturbative theoretical results (solid lines) for the 
two-photon absorption in Na induced by the shaped pulses having a 7r-step spectral phase pattern. 

The results include an integration over the experimental spatial beam profile. The traces show the 
final 4s population 

P4:S ^ 

function of the step position cogtep- The value of P^g is normalized 
by the final population P4s,tl excited by the corresponding transform-limited (TL) pulse. Each 
column corresponds to a different central spectral wavelength Ao=773, 777, and 779.5 nm with a 
different pulse energy in each row (increasing from top to bottom). Each pulse energy corresponds 
to a different transform-limited intensity at the peak of the spatial beam profile _ 
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FIG. 3: Non-pcrturbativc theoretical results (gray thick lines) and 4*'*-order perturbative theoret- 
ical results (black thick lines) for the two-photon absorption in Na induced by the shaped pulses 
having a 7r-step spectral phase pattern. As a 4*'*-order calculations, the perturbative calculations 
include contributions from the 2"*^ and 4*'* orders. The traces show the final 4s population P/^g 
as a function of the step position Wstep- The right-hand y-axis scale is the absolute value of Pa,s- 
The left-hand y-axis scale corresponds to normalized by the final population P45,tl excited by 
the corresponding transform-limited (TL) pulse. Each column corresponds to a different central 
spectral wavelength Ao=773.5, 777, and 780 nm with a different (single-valued) spectral intensity 
7o in each row (increasing from top to bottom). Each /q corresponds to a different transform- 
limited peak intensity Jtl- The weak-field traces are given in the first-row panels [(a)-panels]. For 
comparison, they are also shown in all the other panels (black thin lines) on the TL-normalized 
scale. 
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FIG. 4: Perturbative theoretical results (of 4*'^ order) for the two-photon absorption in Na induced 
by the shaped pulses having a vr-step spectral phase pattern. The traces show the final 4s population 
P/j^s as a function of the step position ujstep- The value of is normalized by the final population 
-f4s,TL excited by the corresponding transform-limited (TL) pulse. Each panel corresponds to 
a different central spectral wavelength Ao=773.5, 777, and 780 nm. The traces given in gray 
thick lines are the real traces (shown also in Fig. [sj at /tl = 2.5 x 10^*^ W/cm^, which is the 
upper intensity limit of the present intermediate-field regime. The traces given in black thick lines 
(without the intermediate- field dip feature at 781.2 nm) are traces that have been calculated at 
/tl = 2.5 X 10^*^ W/cm^ with artificially setting to zero the contribution from the resonance- 
mediated Raman transitions via 7p, i.e., A^^'''^^^ = (see text). For completeness, also presented 
are the real weak-field traces (black thin lines; also shown in Fig. p|. 
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FIG. 5: Theoretical results for different quantities that are introduced by the frequency-domain 
4*'*-order perturbative description. They are discussed in the text for explaining the different 
intermediate- field features. Each panel shows the Raman term A^^°^'^^^^\6) (gray thick lines) and 
the two-photon transition term A^^\lo3s,3s—S) (black thick lines) as a function of the detuning S (see 
text). The detuning values (x-axis scale) are actually given as the normalized values S/Aoo, with 
Auj being the bandwidth of the intensity spectrum of the pulse. The zone of small \6\ around 6 = 0, 
which is the most contributing to yj^Oi)near-res ^^^^ text), is indicated schematically (dashed area). 
Each column corresponds to a different central spectral wavelength Ao=773.5, 777, and 780 nm. 
Each row corresponds to a different pulse shape: (a) transform-limited (TL) pulse, (b) shaped 
pulse with ujstep = ^4^s.3s/'^ (777 nm), and (c) shaped pulse with ujstep = i^o (the corresponding 
central spectral frequency) . For comparison, each panel also shows the A^'^^ {co3s,3s — ^) term of the 
corresponding TL pulse (black thin lines). 
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